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ABSTRACT 

We investigate how the imprint of Faraday rotation on radio spectra can be used 
to determine the geometry of radio sources and the strength and structure of the 
surrounding magnetic fields. We model spectra of Stokes Q and U for frequencies 
between 200 MHz and 10 GHz for Faraday screens with large-scale or small-scale 
magnetic fields external to the source. These sources can be uniform or 2D Gaussians 
on the sky with transverse linear gradients in rotation measure (RM), or cylinders 
or spheroids with an azimuthal magnetic field. At high frequencies the spectra of 
all these models can be approximated by the spectrum of a Gaussian source; this is 
independent of whether the magnetic field is large-scale or small-scale. A sine spectrum 
in polarized flux density is not a unique signature of a volume where synchrotron 
emission and Faraday rotation are mixed. A turbulent Faraday screen with a large 
field coherence length produces a spectrum which is similar to the spectrum of a 
partial coverage model. At low and intermediate frequencies, such a Faraday screen 
produces a significantly higher polarized signal than Burn’s depolarization model, as 
shown by a random walk model of the polarization vectors. We calculate RM spectra 
for four frequency windows. Sources are strongly depolarized at low frequencies, but 
RMs can be determined accurately if the sensitivity of the observations is sufficient. 
Finally, we show that RM spectra can be used to differentiate between turbulent 
foreground models and partial coverage models. 

Key words: magnetic fields - polarization - galaxies: ISM - galaxies: jets - galaxies: 
magnetic fields 


1 INTRODUCTION 

Most Active Galactic Nuclei (AGN) are so distant that we 
can only study their structure using Very Long Baseline 
Interferometry (VLBI). However, the magnetic field in the 
medium surrounding AGN leaves an imprint on the polar¬ 
ized fiux density spectra, allowing us to study AGN using 
polarization measurements. In a magnetized plasma the re¬ 
fractive index for the left-handed and right-handed circu¬ 
lar polarization modes is slightly different, which leads to 
a net rotation of the linear polarization vector between fre¬ 
quencies, also known as Faraday rotation. The amount of 
Faraday rotation of a polarized wave which is detected at a 
wavelength A is equal to Xoba — Xo = RMA^, where Xobs 
and Xo indicate the position angle of the polarized electro¬ 
magnetic wave that is detected and the position angle of the 
wave when it is emitted, respectively. The rotation measure 
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RM is equal to 

^observer 

RM = 0.81 / neB||dl rad m“^, (1) 

J source 

where Ue is the local electron density in units of cm“®, i?|| 
the line-of-sight component of the magnetic field in units 
of plG, and the path length is measured in parsec. If rie, 
i3||, or the path length changes inside or across the source 
then equation[T]applies to each individual region where these 
quantities are the same. In that case the net observed po¬ 
larization vector is a linear combination of the polarization 
vectors which are emitted by volumes with the same RM, 
and the net RM of the source = dXoba/dA'^, which can de¬ 
pend on frequency. 

Depending on the geometry of the magnetic field, its 
strength, coherence length, and the electron density distri¬ 
bution along the line of sight, different parts of the back¬ 
ground source will have different RMs. These regions with 
different RMs leave an imprint on the polarized flux density 
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and position angle spectrum of the source, which can be de¬ 
composed into emission at diff erent RMs using the techniqu e 
of rotation measure synthesis (jBrentiens fc de BruvnlBoOiih . 
In this Paper we investigate how polarization measurements 
can be used to study the properties of synchrotron-emitting 
and Faraday-rotating media in unresolved sources. These 
sources can be distant AGN or structures in nearby galax¬ 
ies which are below the resolution limit of a telescope. 
We develop and analyse polarized flux density and po¬ 
sition angle spectra as a function of frequency, and RM 
spectra, for different geometries of the sources themselves 
and of the magnetic fields in their surroundings. Ana¬ 
lytical models of structured magnetic fields (either large- 
scale oi_small^ale}_Jiave b een p r esented previously by , 
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polarimeters on radio telescopes only provided data for 
very narrow frequency windows. To study the structure 
of magnetic fields in AGN one would typically analyse 
data a t a number of frequ e ncies, as was don e for exam- 
ple by iRossetti etaT _j200^)jjMantoymii_et_^ J 20M ), an d 


more recently by Fames et al.l ~ ^14h . Fames et aLn2014ll . 
Modern radio interferometers like the Jansky Very Large 
Array and the Australia Telescope Gompact Array that 
are equipped with broad-band polarimeters can be turned 
into powerful tools for studying the co mplex structure of 
AGN and th eir associated magnetic fi elds. [Farnsworth et al] 
ll201lll . and lO ’Sullivan et ahl (l2012l l studied compact ra¬ 
dio sources using broad-band data and the t e chniq ue of 
RM synthesis . Building on the w o rk by iBur I ll 19691 and 
[iMbM^ lll99lll , lArshakian fc Be^ ll20li] l and Bernet et al.l 
I 2012tl have modelled polarized flux density spectra for me¬ 
dia that cont ain large-scale and small-sca le magnetic fields. 
Very recently, iHorellou &: Fletcher! ll20l4l presented analyt¬ 
ical models of helical magnetic fields and calculated polar - 
ized flux density spectra from these. IShneider et al.1 (120141 1 
developed analytical models of depolarization by media with 
large-scale or turbulent magnetic fields and used these to in¬ 
terpret their observations of the grand-design spiral galaxy 
M51. 

In our Paper we consider radio sources which are com¬ 
pact compared to the size of the telescope beam, so that 
beam attenuation is not important. Faraday rotation oc¬ 
curs in front of the source, and can be due to either large- 
scale or turbulent magnetic fields. The large-scale magnetic 
fields that we consider produce a transverse linear gradient 
in RM across a uniformly emitting source or a source with 
a Gaussian flux density profile, or they wrap around one 
of the axes of the source. We consider turbulent magnetic 
fields that can be characterised by a single coherence length. 
We develop semi-analytical models for these different source 
types, which we use to calculate polarized flux density and 
position angle spectra as a function of frequency, and RM 
spectra. Then we analyse which features in the frequency 
spectra or RM spectra can be used to identify the geometry 
of the emitting region and the properties of the Faraday- 
rotating screen in front of this region. In our models we 
will consider RMs of up to at least 1000 rad m“^ to reflect 
the range in RMs that have been f ound in AGN cores, see , 
e.g.. lZavala fc Tavlorl ll2003l . l2004ll . lO’Sullivan et al.l ll201lll . 
and iHovatta et aL I 2012ll . Compared to previous studies, we 


consider a wider range of models and develop new models, 
present model spectra for a wider frequency range, include 
depolarization across the finite frequency channels, and anal¬ 
yse these spectra using RM synthesis. The wide frequency 
range that we consider will be accessible with the Square 
Kilometre Array (SKA) which is expected to produce its 
first results around 2020. 

In Section [2] we describe the models that we developed 
that include large-scale or turbulent (small-scale) magnetic 
fields, and we present polarized flux density spectra for these 
source geometries. In Section |3] we calculate RM spectra 
from the Stokes Q and U frequency spectra that we modelled 
for four observing windows. In the Appendix we describe in 
detail how we modelled emitting cylinders and emitting el¬ 
lipsoids with a wrapped-around Faraday-rotating layer, and 
we derive probability density functions of RM for turbulent 
Faraday screens. 


2 MODEL DESCRIPTIONS AND 
FREQUENCY SPECTRA 

When we integrate through the synchrotron-emitting 
medium to calculate the polarized flux density profile of the 
source we make the following assumptions: 

(i) The monochromatic volume emissivity Si, and the 
shape of the emission spectrum are the same throughout 
the source, 

(ii) The emission is synchrotron-thin; synchrotron self¬ 
absorption is not important, 

(iii) The intrinsic position angle of the polarized emission 
is constant throughout the source, 

(iv) Wavelength-independent depolarization is constant 
throughout the source, and 

(v) The bulk velocity of the synchrotron-emitting plasma 
is much smaller than the speed of light. 


The polarization properties of jets w ith relativistic bulk 
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I 201llf . and Forth et al.l (l201ll l. These authors have 
demonstrated that relativistic jets with helical magnetic 
fields inside the source region show complex polarization 
behaviour, for example, the electric vector position angle 
can change direction by 90 degrees between the axis and 
the edge of the jet (‘spine-sheath’ struc tures: see, e.g., 
lAttridee et al.1 [19991 . IPushkarev et ahllioosl l. and the polar¬ 
ized flux density does not have to be symmetric along a 
cross-cut perpendicular to the axis of the jet. In our Paper 
we limit ourselves to sources (which could be jets) with 
non-relativistic bulk velocities to reduce the complexity 
of the models, thereby simplifying our analysis. In all our 
models the Faraday-rotating gas is non-relativistic; for 
Faraday rota tion in a relativisti c gas w e refer the reader to 
the paper bv iBroderick fc LoebI (l2009ll . 

We will use the ‘intrinsic polarization fraction’ po to re¬ 
fer to the polarization fraction that is measured by a hypo¬ 
thetical observer between the source of the emission and the 
Faraday-rotating foreground screen, while we ourselves mea¬ 
sure the polarized emission only after it has passed through 
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the foreground screen. This intrinsic polarization fraction 
includes wavelength-independent depolarization effects that 
are not caused by Faraday rotation. Dividing the surface- 
integrated monochromatic polarization vector by po times 
the Stokes I spectrum removes spectral index effects. We 
choose the orientation of the Stokes Q, U coordinate sys¬ 
tem such that the source emits only in Stokes Q (before the 
emission passes through the Faraday-rotating foreground 
medium). 

The observed monochromatic polarization vector at fre¬ 
quency V is given by 


T*obs 


f Pern dxdy(X&) 

J source 

/ oo 

Pem (RM) dRM, (3) 

- OO 


where Pem {x, y, z) indicates the polarization vector that is 
emitted at position (a:, y, z) within the source, Pem (RM) 
indicates the polarization vector that is emitted at a single 
RM, and ‘c’ the speed of light. Throughout this Paper we 
will use a Cartesian coordinate system x, y, z where x points 
towards the observer, and y and z lie in the plane of the sky 
such that 2 : points along the major axis of the source and y 
along the minor axis. The major axis of the source can be 
inclined by an angle 9 with respect to the plane of the sky; 
in that case y and 2 point along the projection of the minor 
and the major axis of the source on the sky, respectively. 

We will use analytical expressions or numerical mod¬ 
elling to determine Pobs p)- Then we calculate the net po¬ 
larization vector of each frequency channel by integrating 
over the response function of that channel. For channel j 

POO 

Pobs(^'j) = / Wj (u) Pahs{l^)du . (4) 


The response function of channel j, Wj (v), expresses how 
sensitive the observations are to each frequency in the ob¬ 
serving band. The integral of Wj {v) out to ± infinity is equal 
to one. We will simulate frequency channels with a top-hat 
response in frequency, which is equal to 1 15v inside the chan¬ 
nel and zero outside it, where Sv = — Vj is the channel 

width. We will refer to the channel-averaged polarization 
vectors as Pobs {v\Sv). The polarized flux density that the 
background source emits is independent of the structure of 
the Faraday screen, therefore 

/ OO 

|Pem(RM)|dRM = constant. (5) 

-OO 


In our simulations the background source emits 1000 units of 
polarized flux density. The physical interpretation of equa¬ 
tion [ 5 ] is that structure in the Faraday screen redistributes 
polarized flux density over different RMs. As a result sources 
that emit over a wider range of RMs have a lower peak po¬ 
larized flux density in the RM spectrum. 


2.1 A foreground medium with a large-scale 
magnetic field 

In this section we simulate transverse linear RM gradients 
in front of a source that emits uniformly across its surface 
(§ 1 1^.1. Ill or has a Gaussian emission profile (§ 12.1.211 . Such 
RM gradients can occur when the source of the emission 
is embedded in a large ionized halo. If the source is small 


compared to the radius of the halo then any inclined jet 
will show a transverse linear RM gradient across its surface 
because of the increase in path length through the halo. 
Next we model sources where the geometry of the Faraday- 
rotating foreground and the geometry of the emitting region 
are closely connected. In S I2.1.3l we consider a magnetic field 
that wraps around a cylinder which emits polarized radio 
waves; the magnetic field inside the emitting cylinder points 
along the major axis of the cylinder. In S I2.1.4l we generalize 
the shape of the emitter to an ellipsoid that can have any 
axis ratio and any orientation with respect to the line of 
sight. Azimuthal m agnetic fields have been o bserved in the 
radio jets of 3C273 (lAsada et al.ll20o3j. M87 (lAlgaba et al. 
201 s ) ■ and in other radio jets fe.g.. iReichstein fc Gabuzda 
2012 , iGabuzda et al.l[2014l ~l. even though in those cases the 
synchrotron-emitting plasma moves down the jet at rela¬ 
tivistic speeds, while for the cylinders and ellipsoids that we 
model we consider plasmas with non-relativistic bulk veloc¬ 
ities. The expressions we derive are valid for a single source; 
if a source consists of multiple components then the result¬ 
ing RM spectrum is simply the vector sum of the complex 
RM spectra of the individual components. 

To avoid unnecessary calculations, when we use numer¬ 
ical integration we increase the (fractional) numerical ac¬ 
curacy of the integral at the highest frequencies, where the 
integrated polarized flux densities are highest. This guaran¬ 
tees that the integrated flux densities at all frequencies are 
accurate down to at least 10“® flux density units, and that 
we do not spend too much time numerically integrating at 
the lowest frequencies. We start our simulation always at 
the highest frequency, where depolarization is minimal, and 
for this frequency we set the fractional numerical accuracy 
e = 10“^°. Then we move to lower frequencies, and we let e 
for channel j increase as 


G+i 


]^Q-floor[logio(abs(f>o|,s,j)/OOl)-|-5] 


( 6 ) 


up to 10“®. When numerically integrating non-monotonic 
functions we checked that the numerical integral was calcu¬ 
lated with a sufficiently small (accurate) e. In ? l2.1.4l we will 
use a different method to fix the numerical accuracy with 
which polarization vectors are calculated for each of the fre¬ 
quency channels. The accuracy of our numerical method is 
sufficient to calculate RM spectra in Section |3] down to at 
least 0.01 flux density units (-50 dB). 


2.1.1 Uniform source with a transverse linear RM 
gradient 

The net monochromatic polarization vector Pobs (v') for a 
uniform source of emission inside which RM increases or 
decreases linea rly wi t h phy sical depth has been calculated 
analytically bv iBurnI (1 196611 : 

Pohs{v) = (po X f) sine (ARM (c/iz)^) , (7) 

where ARM = RMmax — RMmin is the range in RM over 
which the source emits, and RMc = 0.5RMmin + 0.5RMmax 
is the mean RM of the emission. We define sine ( 2 ;) = 
sin (a:)/a:. Such a source in which emission and Faraday 
rotation are mixed has become known as a ‘Burn slab’. 
Equation [ 7 ] also describes the spectrum of a uniformly emit¬ 
ting source which has a linear RM gradient in front of it 
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(E.g.. [Sokoloff et n]|l998l . and IZavala fc Tavlorj|2004) . Fig.[T] 
shows the absolute value of the simulated polarized flux den¬ 
sity Poha (v'jSu), scaled with the assumed polarized frac¬ 
tion times the flux density spectrum in Stokes I, of uniform 
sources that lie behind linear gradients in RM. We inte¬ 
grated Pobs (iz) across frequency channels using Romberg’s 
method. 

Equation [7] expresses how the emission at different RMs 
leads to depolarization (the sine term) while the polarization 
vectors at different frequencies show Faraday rotation with 
a net RM equal to the mean RM of the emission. If RMc 
= 0 rad m“^, then the observed polarization vectors do not 
show net Faraday rotation, and all polarized emission will 
be in Stokes Q, Stokes U will be 0. 

The condition that |Pem (RM)| is the same for all RM is 
much more general than one might think at hrst. There are 
no constraints on the geometry of the source, nor on the in¬ 
clination of the source with respect to the plane of the sky. 
The source can even consist of multiple components. The 
RMs across the source do not have to increase or decrease 
monotonically, and the RM gradient can have any orienta¬ 
tion with respect to the major axis of the background source. 
For example, a cylindrical source with a linear RM gradient 
along its major axis is described by equation [7] Even the 
spectrum of a source with an azimuthal magnetic field can 
be approximated by equation [7] if the emission prohle along 
the RM gradient of the source is approximately constant. 
Unfortunately, because all these geometries have the same 
frequency and RM spectra it is not possible to determine 
the source geometry without additional information. 


2 . 1.2 


Leah^et^ 


Sokoloff et al.l 


Gaussian source with a transverse linear RM 
gradient 

I Johnson et ahl 
calculated 


(H), 

monochromatic 


and 

polar¬ 


ization vectors for a uniformly emitting source with a linear 
gradient in RM in front of it, which is observed with a Gaus¬ 
sian beam. We note that by re-ordering the factors inside 
the integrals that these authors solved the same solution 
is found for a source with a 2D Gaussian profile that lies 
behind a Faraday screen with a transverse linear gradient 
in RM. The Gaussian source measures Oy x cr^ standard 
deviations along its minor and major axis, respectively, and 
the RM gradient is described by the function RM(t/, z). The 
observed polarization v ector at a single fre quency is then 
given by equation 40 in ISokoloff et al.l (Il998ll : 

(Po X I) 


/ oo rexi 

/ ® 

-oo J — oo 

e2®MG.D(c/G"dt/d2 




= (po X 7)e 


- 2ARM’^(c/i^)'‘ + 2iRMc(c/i^)^ 


( 8 ) 


where 


ARM^ = 


/9RM V /9RM 
' -0-„ I -I- I -(Tj 


V 9y 


V dz 


is the total change in RM across the Gaussian source, and 
‘RMc’ is the RM at the centre of the Gaussian source on 
the sky, at coordinates {yo,zo). ARM^ can be expressed in 
terms of the RM difference across the FWHM of the major 
and minor axes of the source by using a = FWHM / V 81n2. 
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Figure 1. Broad-band polarized flux density spectra between 
200 MHz and 10 GHz, scaled with the intrinsic polarization 
fraction times the Stokes I spectrum, for Gaussian background 
sources that illuminate a linear gradient in RM. The source emits 
1000 flux density units. The RM differences across the FWHM of 
the major axis of the Gaussian are 10, 250, and 1000 rad m~^ (red 
to blue). The grey lines in the background show spectra for uni¬ 
form sources from ^ I2.1.1l with a half-width equal to three times 
the standard deviation of each of the simulated Gaussian sources. 
These uniform sources emit over RM ranges of 13, 318, and 1274 
rad m“^. The spectrum of a Gaussian source and the matching 
spectrum of a uniform source intersect at approximately 300 flux 
density units. 


We numerically integrate equation [8] over the widths 
of the individual frequency channels using Romberg inte¬ 
gration. Because \Pohs{i')\ decreases monotonically with 
decreasing u, and because the channel width Su is con¬ 
stant, \Pohs decreases monotonically with decreas¬ 

ing u. Starting at the highest observing frequencies, we 
stopped integrating equation [8] at the frequency where 
I-Pobs Su) \ / (po X I) was < 10~^, which is below the limit 
we consider in Section [3] for calculating RM spectra. 

Fig. □ shows frequency spectra of Gaussian sources 
with transverse linear RM gradients between 10 and 1000 
rad m~^ PWEM"*^. Gaussian emitters with even a small RM 
difference across the source (> 125 rad FWHM“^) are 
almost completely depolarized at frequencies below 1 GHz. 
This is not the result of depolarization across the frequency 
channels, which we checked by simulating observations with 
frequency channels of 0.1 MHz. When the RM gradient is 
small, the polarized flux density spectrum shows a sharp 
drop over a narrow frequency range. For such sources it is 
vital to match the frequency window of the observations to 
the RM gradient one is looking for. If the observing fre¬ 
quency is too low the source will be severely depolarized 
and therefore difficult to detect. On the other hand, if the 
observing frequency is too high the source will not show 
any depolarization, and will be indistinguishable from other 
sources with small RM gradients which also do not show 
depolarization. 

If a source has a Gaussian profile along its minor axis 
and a uniform profile along its major axis, and illuminates a 
linear RM gradient, then the orientation of the RM gradient 
determines if the spectrum we observe follows the sine pro¬ 
file from § 12.1.11 (RM gradient along the major axis), that 
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of a Gaussian source (RM gradient along the minor axis), 
or a combination of the two. If a source consists of uniform 
and Gaussian components that from our perspective all lie 
behind the same linear gradient in RM then the frequency 
and RM spectra are a combination of the spectra which we 
modelled in § 12.1.11 and this section. At high frequencies 
the uniform source depolarizes at roughly the same rate as 
a Gaussian source of the same extent, i.e. the long axis of 
the uniform source has a half-length equal to approximately 
three times the standard deviation of the Gaussian source. 
The grey lines in the background of Fig. [T] show the spec¬ 
tra of uniform sources with the same extent as the Gaus¬ 
sian sources that we modelled; the lines for each matching 
pair of spectra intersect at about 300 flux density units. At 
high frequencies the difference between the spectra of the 
two source types becomes more pronounced when the RM 
gradient across the Gaussian and uniform source becomes 
steeper, while at low frequencies the re-brightening of spec¬ 
tra is a tell-tale sign for uniform sources, which is not found 
in Gaussian sources. If the uniform source is larger than the 
matching Gaussian source it will be more depolarized than 
the Gaussian source, and vice versa. 

2.1.3 Cylindrical source with an azimuthal magnetic field 

For this geometry, Faraday rotation occurs in a boundary 
layer between an inner, emitting cylinder of radius R and 
an outer, coaxial, cylinder of radius R'. The cylinders can 
have any inclination with respect to the plane of the sky. In 
Appendix m we explain in detail how we model frequency 
spectra for this source type; there we also show that the 
frequency spectra that we model are independent of the in¬ 
clination angle. This implies that the inclination angle of the 
cylinder can not be derived from the frequency or RM spec¬ 
trum of a source. For this simple geometry the only input 
parameters are the maximum RM of the source, RMmax, 
which is found at the edges of the emitting cylinder, and 
the intrinsic polarized flux density that is emitted by the 
source, which we normalize to 1000 units of polarized flux 
density. Our model of the cylinder also acts as a test case for 
the more complex model of an emitting ellipsoid that we will 
discuss in S I2.1.4I We found that a few of the 9800 frequency 
channels that we simulated do not decrease monotonically in 
polarized flux density with decreasing frequency; this could 
be the result of our numerical integration scheme. Because 
the difference in polarized flux density is at most 0.2 units of 
polarized flux density, and often at least one order of magni¬ 
tude smaller, and because only very few frequency channels 
are affected, the RM spectra that we calculate from our sim¬ 
ulation are accurate down to at least 0.01 units of polarized 
flux density. 

Fig. [2] shows polarized flux density spectra for different 
values of RMmax, if Faraday rotation occurs in a boundary 
layer with a thickness equal to 1% or 10% of the radius of the 
emitting cylinder. Increasing the thickness of the boundary 
layer leads to more depolarization, which has the effect of 
spreading out the curves that correspond to different RMmax 
in Fig. [21 This makes it easier to tell the curves for the 
different RMmax apart, and it makes it easier to determine 
RMmax more accurately from observations that cover a wide 
frequency range. Fig. [2]also shows that spectra of a cylinder 
with a thick Faraday-rotating layer and a small RMmax can 
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Figure 2. Frequency spectra for a cylinder of polarized emission 
with an azimuthal magnetic field which wraps around the major 
axis of the cylinder. The solid lines indicate cylinders with RMmax 
of 25, 500, and 1500 rad m~^ (red to blue) and a Faraday-rotating 
layer with a thickness equal to 10% the radius of the emitting 
cylinder; the dashed lines show results for the same RMmax but 
a Faraday-rotating layer with only 1% the radius of the emitting 
cylinder. The grey spectra in the background are for Gaussian 
sources with a transverse linear RM gradient of 7, 145, and 436 
rad m“^ FWHM“^ across its major axis. 

be very similar to spectra of a cylinder with a thin boundary 
layer and a large RMmax, which is illustrated by the blue 
dashed curve which overlaps with the green solid curve. This 
will make it difficult to tell the thickness of the boundary 
layer from the observed frequency spectrum. 

In Fig. (21 we also compare the frequency spectra for a 
cylinder with a thick Faraday-rotating layer with frequency 
spectra for a Gaussian source with a transverse linear RM 
gradient across its major axis, which we modelled in § 12 .1.21 
We chose the values for the RM gradients in front of the 
Gaussian source such that frequency spectra of the cylin¬ 
der and the Gaussian source intersect at 500 polarized flux 
density units. The comparison between the coloured and 
greyscale frequency spectra shows that it could be possible 
to tell the two source types apart by accurately measuring 
the difference in curvature of the frequency spectra, which 
becomes easier if the amount of Faraday rotation across the 
Gaussian source or across the minor axis of the cylinder is 
large. 

2 . 1.4 Ellipsoidal source with an azimuthal magnetic field 

In this model, Faraday rotation occurs between an inner 
and outer ellipsoid, where the magnetic field in the Faraday- 
rotating layer wraps around the major axis {z axis) of the in¬ 
ner ellipsoid. These ellipsoids have a common coordinate sys¬ 
tem that we introduced at the beginning of this Section, but 
their extents can be chosen independently. In Appendix [B] 
we describe in detail how we modelled this source type. The 
ellipsoid is described by many parameters; we reduce the di¬ 
mensionality of the parameter space by considering the case 
of a spheroid, which has a circular cross-section of radius R 
perpendicular to the major axis 2 . The spheroid can be de¬ 
scribed by specifying the axis ratio C/R between the length 
of polar ( 2 ) axis C and its radius at the equator R, the in- 
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Figure 3. Stokes Q profiles (in arbitrary flux density units) at 
1.4 GHz when the surface of a spheroid with axis ratio 1:1;2 and 
a thick Faraday screen is projected on the sky. Different rows 
indicate different inclination angles d of the major axis of the 
spheroid, while different colums show configurations with RMmax 
of 20, 200, and 1000 rad m“^, respectively. 


clination 9 of the spheroid with respect to the plane of the 
sky, the thickness of the Faraday-rotating layer, and RMniax- 
Fig. [3] shows the pattern of Stokes Q across the surface of a 
spheroid which has different inclinations with respect to the 
plane of the sky and different RMmax- 

Fig.Ushows frequency spectra for spheroids with differ¬ 
ent axis ratios, inclination angles, RMmax, and thicknesses of 
the layer of Faraday-rotating layer. The numerical accuracy 
of these spectra is < 10 units of polarized flux density; in Ap- 
pendixl^we explain how we derived this value. If the source 
lies in the plane of the sky the different axis ratios that we 
tested produce the same frequency spectrum to within a few 
units of polarized flux density at 500 MHz; therefore we only 
show one of the three panels for this inclination angle. Sim¬ 
ilar to what we noticed for the cylindrical source, increasing 
the thickness of the Faraday-rotating layer spreads out the 
frequency spectra. Changing the inclination of the spheroid 
clearly has a larger impact on the spectra shown in Fig. 0] 
than changing the axis ratio of the spheroid. Therefore it 
is easier to determine the inclination of the spheroid than 
its axis ratio. If the source is inclined with respect to the 
plane of the sky frequency spectra can show re-brightening 
at low frequencies, with secondary maxima which reach the 
same height as the secondary maxima of the uniform source 
discussed in § 12.1.11 As we concluded in § 12.1.11 this oc¬ 
curs when the emitted polarized flux density is the same for 
all RM at which the source emits. Needle-like sources with 
large C/R axis ratios do not show a strong amplitude differ¬ 
ence between the secondary and higher-order maxima and 
minima. 

At high frequencies the spectrum of a spheroid can 
be approximated by the spectrum of a Gaussian source 
which we modelled in S I2.1.2I Accurate measurements at low 
and intermediate frequencies can detect the secondary and 
higher-order maxima which are present in the spectrum of 
an inclined spheroid but not in the spectrum of a Gaussian 


source; such measurements might even distinguish between 
an inclined spheroid and a uniform source with a transverse 
linear RM gradient (S I2.1.1II . 

The spheroids that we modelled do not show a change 
in position angle with frequency beyond the numerical accu¬ 
racy of the models. Therefore, the net RM of the spheroids 
is zero, which we also found for the models with large-scale 
magnetic fields which we considered previously. 

2.2 A turbulent foreground medium 

Finally, we consider depolarization due to a turbulent mag¬ 
netic field in front of a uniformly emitting source. We inves¬ 
tigate how the coherence length and magnetic field strength 
affect the polarized flux density spectrum as a function of 
frequency. 

Depolarization by a turbulent foregr o und h as been in¬ 
vestigated by a number of authors. iBurnI (Il966l) considered 
a foreground screen with RMs that are drawn from a Gaus¬ 
sian distribution with variance All cells in this model 

have the same volume equal to the c ube o f the coherence 
length of the magnetic field. iTribbld lll99ll) expanded this 
analysis to in c lude s tructure in RM on a spectrum of scales. 
iMurgia et al.l ll2004l ) modelled the polarization properties of 
galaxy clusters in 3D using a power-law to describe structure 
in the magnetic field, and showed that the RM dispersion as 
a function of impact parameter, ctrm (R±), can be approx¬ 
imated well by a model wit h a single charact eristic length 
scale of the magnetic fielcfl iFanti et al.l (|2004| ) investigated 
how in Compact Steep Spectrum sources the RM dispersion 
follows the King profile ctrm (Rx) oc (l -I- 
iH in 3 119621: also known as a /3 profile) w here Rc is the 
core r adius of the galaxy under investigation. I^ssetti et al.l 
ll2008l) proposed a model where the turbulent foreground 
only covers a fraction of the emitter in the background. In 
their model, the observer sees a combination of the Burn 
depolarization model and a contribution by the background 
source that is not covered by the turbulent foreground. The 
RMs in their model are drawn from a Gaussian distribution, 
and there are at most two Fa r aday- rotating screens with tur¬ 
bulent layers. iHovatta et al.l (l2012l) modelled depolarization 
by a turbulent Faraday screen with a Gaussian probability 
density function (pdf) of RM in front of a radio jet. They 
concluded that the depolarization of isolated jet components 
in their sample of AGN can be explained by a small number 
of sightlines passing through such a Faraday screen. 

We model a circular source with a diameter of 25 pc that 
emits 1000 units of polarized flux density uniformly across 
its surface. The results we derive apply also to sources which 
are not circular, as long as the surface area of those sources 
is the same as the surface area of the circular source we con¬ 
sider. The coherence length of the magnetic field is allowed 
to vary continuously; therefore the number of turbulent cells 
does not have to be an integer. The intrinsic polarization 
angle of the emission is set to zero degrees throughout the 
emitter. We model a Faraday-rotating screen in front of this 
source that consists of one layer of turbulent cells (§ 12 . 2. Ill 

^ This length scale is equal to J)j°° ui|| (R)dA/ui|| (0), where 1 U|| (i?) 
is the spherically-averaged autocorrelation function of the mag¬ 
netic field component along the line of sight. 
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Figure 4. Frequency spectra for different axis ratios of the spheroid (columns), inclination angles 9 (rows), and RMmax of 25, 500, and 
2500 rad (red to blue). One in every 25 channels is shown. The axes of the outer spheroid are scaled versions of the axes of the 

inner spheroid, where the scale factor is either 110 per cent, producing a thick layer of Faraday-rotating material (solid lines) or only 101 
per cent, producing a thin layer (dashed lines). The axis ratios of each column are indicated above the top panel in each column. The 
three panels for a spheroid which lies in the plane of the sky {6= 0°) are almost identical; we show only one of the panels (see the text 
for details). The grey lines in the background show frequency spectra for a Gaussian source with a transverse linear RM gradient, which 
we modelled in § 12.1.21 The RM gradients in front of these Gaussian sources were chosen such that the spectrum of a spheroid and the 
spectrum of the matching Gaussian source intersect at 500 flux density units. 


or multiple layers of cells (§ 12.2.211 . The electron density tie 
is constant throughout the screen, and equal to 10 cm“®. 
We investigate magnetic field strengths Bturb of 1, 5, 10, 
25, and 50 /rG, while the direction of the magnetic field 
is drawn independently between turbulent cells. One can 
choose a different value for the electron density or magnetic 
field strength; as long as the product of the electron den¬ 
sity times the magnetic field strength is the same as in our 
simulation one will find the same frequency spectra and RM 
spectra. 


2.2.1 A single layer of cells 

For a single layer of turbulent cells the observed monochro¬ 
matic polarization vector is equal to 

.0. S. / /. \ 

PobPv) = ^ ^y^ Pem(i/,2)d2/d2j (9) 


‘l.o.s.’ is shorthand for ‘line of sight’, there are A^i.o.s. sight¬ 
lines, and the coordinate axes y and 2 point along the 
minor respectively the major axis of the emission region. 
Because the source that we model emits uniformly across 
its surface, the integral in equation [9] can be replaced by 
P/ (po X /) = 1000 flux density units divided by the num¬ 
ber of lines of sight through the Faraday screen. This sim¬ 
plifies equation |9] to 


P obs (n) 


P/(poX I) 

Nl.o.s. 



1=1 




( 10 ) 


In Appendix o we show that if the magnetic field direc¬ 
tion is drawn randomly, such that every point on the unit 
sphere has equal probability of being drawn, then pdf(RM) 
is uniform between ± the largest possible RM of the screen, 
RMmax- In our simulation the maximum RM of a single tur¬ 
bulent cell is given by 


RMmax = 


O.SlBturbUeyTt (25/2 pc)^ /All.o.s. - 


( 11 ) 
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Figure 5. The equivalent of Fig.[T]for a single layer of 10 turbulent cells (left panel) with RMs between ± 569 rad m“^, and a single layer 
with 324 turbulent cells (right panel), with RMs between ± 100 rad m“^. The greytones indicate the different confidence intervals, which 
contain 95% (dark grey), 75% (white) and 50% (black) of the simulated ensemble. These RMs are drawn from a uniform probability 
density function for RM, and we repeated this process 2000 times to calculate the confidence intervals. The amplitude of the turbulent 
magnetic field is 10 /rG, the free electron density 10 cm“®, and only one in 25 of the 1-MHz channels is shown. The grey dashed line 
indicates how a single frequency channel is depolarized by the largest possible RM in the Faraday screen, RMm ax. The red line shows the 
prediction for the length of the monochromatic polarization vector by the depolarization model bv lBuriJ lll96fih who assumed a Gaussian 
pdf(RM). We calculated the standard deviation of the Gaussian analytically from the pdf(RM) of the monolayer. The red dashed line 
indicates the root-mean-square monochromatic polarized flux density level that is expected for a two-dimensional random walk model. 



2000 4000 6000 8000 10000 

Frequency [MHz] 



Figure 6. The equivalent of Fig. |5] for a turbulent Faraday screen which consists of eight layers of cells. In the panel on the left 
^turb = 25 fiG divided by 8 (the number of layers), which produces RMs between ± 1422 rad In the panel on the right each 

layer in the Faraday screen contains 2022 cells; in this case each cell has a magnetic field strength of 25 fiG. 


To simulate a Faraday screen with a single layer of turbu¬ 
lent cells we calculate the monochromatic polarized flux den¬ 
sity (equation 1101) by drawing one RM per sightline from a 
pdf(RM) that is uniform between ± RMmax. We repeat this 
process 2000 times to build up an ensemble of Monte Carlo 
realisations. The frequency channels that we simulate are 
wide enough that position angles can change by more than 
7 T radians across a single channel. To accurately calculate the 
net polarization vector of discrete channels we calculate for 
each channel how much the position angle changes due to 
the turbulent cell with the largest (absolute) RM. We evalu¬ 
ate equation [To] at forty regularly-spaced frequency intervals 
for each 27t revolution in position angle, and we apply the 


trapezium rule to each of these intervals to calculate the net 
polarization vector across each channel. 

Fig. [5] shows polarized flux density spectra for a turbu¬ 
lent monolayer with 10 cells across the surface of the emitter 
(i.e. a field coherence length of 7 parsec) in the panel on the 
left, and 324 cells (a coherence length of 1 parsec) in the 
panel on the right. The greytones indicate the most com¬ 
pact confidence intervals that contain 50%, 75%, and 95% 
of the Monte Carlo realisations (coloured with dark grey, 
white, and light grey, respectively). As the grey dashed line 
in each panel indicates, depolarization across a single chan¬ 
nel becomes important only at frequencies below ~ 1 GHz. 

The panel on the left in Fig. [5] shows that if only a 
few sightlines pass through the turbulent foreground then 
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one will measure a significant polarized flux density even at 
low frequencies. The scatter in the measured polarized flux 
densities is considerable: 50% of the realisations lie further 
than ~ 75 flux density units from the median value. If many 
sightlines pass through the turbulent foreground then the 
resulting frequency spectrum looks similar to that of a uni¬ 
form source with a transverse linear gradient in RM fS I2.1.1ll 
if both sources emit over the same range in RM (panel on 
the right of Fig. This happens because both source types 
have the same pdf(RM), which is sampled discretely by the 
turbulent foreground and continuously by the linear RM gra¬ 
dient. Increasing RMmax leads to strong depolarization even 
at high frequencies. If both the field strength and the num¬ 
ber of sightlines through the Faraday screen are increased, 
in such a way that RMmax is kept fixed at 100 rad m“^, 
then pdf(RM) is sampled more continuously, and the sinc- 
like spectrum in the panel on the right of Fig. [5] becomes 
smoother. 

The solid red line in both panels of Fig. [5] indicat es the 
predic tion for a Gaussian pdf(RM), as modelled by iBurnI 
lll96(i i. for which we calculated the standard deviation an¬ 
alytically from the uniform pdf(RM) that we simulated. At 
high frequencies the resemblance between the spectra pro¬ 
duced by the Burn model and by the turbulent monolayer 
that we simulated is remarkable. However, at low and in¬ 
termediate frequencies our model predicts a much higher 
polarized flux density than the Burn model. 

To calculate the monochromatic polarized flux density 
the Burn model assumes a continuously sampled pdf(RM), 
which requires a very large number of sightlines through 
the turbulent Faraday screen. The polarization vectors from 
these sightlines align at the highest frequencies that we sim¬ 
ulated, and we detect all 1000 units of polarized flux density 
which the source emits. At lower frequencies Faraday rota¬ 
tion becomes important, the different polarization vectors 
become misaligned because of their different RMs, and this 
leads to depolarization. Because the Burn model assumes 
that the number of sightlines through the turbulent screen is 
very high, depolarization is complete. However, in our model 
of the turbulent screen there are far fewer sightlines than the 
Burn model assumes. All turbulent cells produce polariza¬ 
tion vectors of the same length because we assumed that 
the source emits uniformly across its surface, and because 
depolarization by a large RM across individual frequency 
channels only becomes important at very low frequencies, 
as indicated by the dashed lines in Fig. [S] If these polar¬ 
ization vectors would furthermore have random orientations 
then their behaviour can be described as a two-dimensional 
random walk. The root-mean-square (rms) length of the net 
(summed) polarization vector is equal to the length of the 
individual polarization vectors, |Pem| /(po x f^)/A^i.o.s., times 
pNi.o.s., which describes the rms distance from the origin in 
a two-dimensional random walk using unit vectors. We indi¬ 
cated this rms length of the net monochromatic polarization 
vector with a red dashed line in Fig. [5l It is clear that at 
low and intermediate frequencies a random walk model with 
a small number of sightlines produces a much larger (rms) 
polarized flux density than the Burn model. The red dashed 
line lies at the high-end of the black confidence interval in 
all panels. This, combined with the fact that the simulated 
spectrum shows structure while the spectrum for a random 
walk model does not, indicates that the random walk model 


should only be used for obtaining a rough estimate of the 
expected polarized flux density level at low and intermediate 
frequencies. 

2.2.2 Multiple layers of cells 

While the RM distribution of a single layer of turbulent 
cells can be modelled easily using a Monte Carlo simula¬ 
tion, drawing RMs for each of the cells individually becomes 
increasingly more computationally expensive if a Faraday 
screen consists of multiple layers. Instead, in Appendix [C] 
we determine pdf(RM) analytically for up to twenty layers 
of turbulent cells, and we draw W.o.a. RMs from this pdf, one 
RM for each sightline. We draw RMs from the appropriate 
pdf(RM) using the tech nique of rejection s ampling, which is 
described in section 7 in I Press et"^ lll992l l. The monochro¬ 
matic polarization vector Fobs (v) can then be calculated 
from equation [la using the RMs we draw from pdf(RM) 
as RMj. We use the same method as for the monolayer to 
integrate over the finite width of the frequency channels. We 
assume that different layers are not shifted with respect to 
each other, so that cells in different layers lie on top of each 
other. 

Fig.|6]shows frequency spectra for Faraday screens with 
eight layers of turbulent cells. Spectra for turbulent screens 
consisting of one, two, four, and eight layers of turbulent 
cells behave in a similar way: there is a gradual drop-off at 
the highest frequencies, followed by a flattening at the low¬ 
est frequencies. Spectra of screens that consist of one or two 
layers of turbulent cells show secondary and higher-order 
maxima (Fig. [SJ; such features are absent from turbulent 
screen which consist of at least four layers of cells. At high 
frequencies the spectrum of a Faraday screen which con¬ 
sists of many layers of turbulent cells closely resembles the 
spectrum of a screen with a Gaussian pdf(RM), shown as a 
red solid line in Figs [S] and This similarity of the spectra 
is particularly striking if the number of sightlines through 
the Faraday screen is large, as illustrated by the panel on 
the right of Fig. Surprisingly, pdf(RM) which are clearly 
not Gaussian (Appendix [nj can produce frequency spectra 
which are similar to the spectrum of a Gaussian pdf(RM). 

The shape of pdf(RM) is not important to describe a 
random walk process; what matters is that the direction of 
each polarization vector is drawn independently from the 
same parent distribution as the other vectors. Therefore we 
can apply the random walk model that we introduced in 
§[123] also to turbulent Faraday screens which consist of 
more than one layer of cells. The red dashed line in Fig. O 
indicates the rms length of the net monochromatic polariza¬ 
tion vector. As we found for the monolayer, if the number of 
sightlines through the Faraday screen is small our random 
walk model predicts a higher polarized flux density at low 
and intermediate frequencies than the Burn depolarization 
model. The predicted polarized flux density lies at the high- 
end of the black confidence interval, and should be used as 
a rough estimate of the actual polarized flux density. 

The spectra of many types of astrophysical sources can 
be described at low and intermediate frequencies by the ran¬ 
dom walk model we propose. Models where a background 
source is only partially covered by a thick layer of Faraday- 
rotating turbulent cells have been used in the past to ex¬ 
plain the depolarization behaviour of compact AGN (e.g.. 
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iRossetti et al.ll2008f) and of source s with Mg II absorbers in 
the foreground | Bernet et al ][2 m 3) at low radio frequencies. 
The amplitude of the monochromatic polarized flux density 
in these models is given by 


Table 1. FWHM of the RM spread function for each of 
the frequency windows for which we calculate RM spectra. 
These FWHM values were calculated using equation 61 in 
iBrentiens fc de BruvnI ll2005h . 


\Poh3 {v)\ / (po X /) = /cexp (-2ctrm (c/j/)'^) + (1 - /c)(12) 

where fc indicates the fraction of the polarized flux den¬ 
sity which passes through a Faraday screen with a Gaus¬ 
sian pdf(RM); this Gaussian pdf(RM) has a standard de¬ 
viation (Trm. For a uniformly emitting source fc is equal 
to the surface area of the background s ource that is cov¬ 
ered by this Faraday-rotating screen. Also iBur ^ lll966ll con¬ 
sidered the possibility that there are only a few Faraday- 
rotating ‘clouds’ along the line of sight towards the back¬ 
ground source, so that the background source is only par¬ 
tially covered by the Faraday screen. Partial coverage models 
predict a slow drop-off when going from high to intermediate 
frequencies, which levels off towards even lower frequencies; 
our Monte Carlo model shows a similar behaviour in Figs.[S] 
and[6l The polarized flux density of the plateau in the spec¬ 
trum of a partial coverage model can be predicted also by a 
2D random walk model of the polarization vectors. However, 
partial coverage models do not predict the sharp drop-off in 
polarized flux density that can be seen in our models at 
very low frequencies, which is the result of depolarization 
across individual frequency channels. Observations at these 
frequencies can distinguish between a partial coverage model 
and the Monte Carlo model we proposed. In S I3.2.2l we show 
that RM spectra can also be used to tell these models apart. 


2.3 Comparison between models with large-scale 
and turbulent magnetic fields 

In § 12.11 we showed that the spectra of all the models with 
large-scale magnetic fields which we considered can be ap¬ 
proximated at high frequencies by the spectrum of a Gaus¬ 
sian source with a transverse linear gradient in RM. Spec¬ 
tra of sources with turbulent Faraday screens can be ap¬ 
proximated at high frequencies by the spectrum of a tur¬ 
bulent screen with a Gaussian pdf(RM) which is pierced by 
many sightlines (§ mi)- In fact, such a turbulent Faraday 
screen produces the same spectrum as the Gaussian source 
from §[131 The Gaussian pdf(RM) of this turbulent Faraday 
screen is sampled continuously, therefore 

2 

(po X 7)pdf (RM)e^®“(‘’/") dRM 

- OO 

\p2fn.a J - OO 

= (13) 

as shown bv iBur 3 l|l96fil. equation 21 in his paper). Equa- 
tion[T3]depends on frequency in the same way as equation^ 
modulo a term in the latter equation which expresses Fara¬ 
day rotation of the emission. 

Determining from observations at high frequencies 
whether a source has a large-scale or turbulent magnetic 
field is therefore not straightforward. One can tell the dif¬ 
ferent models apart by accurately measuring the shape of 
the spectrum at high frequencies, and by including mea¬ 
surements at intermediate and low frequencies if available. 


Frequency band 
(MHz) 

FWHM 
(rad m~^) 

350 - 900 

6 

950 - 1760 

54 

1300 - 3100 

87 

5000 - 7000 

2158 


One advantage of the similarity of spectra at high frequen¬ 
cies is that the model of a Gaussian source can be used to 
calculate the approximate shape of the spectrum of a more 
complex source very quickly. 


3 ROTATION MEASURE SPECTRA 

In this Section we simulate RM spectra for four frequency 
bands: 350-900, 950-1760, 1300-3100, and 5000-7000 MHz, 
to understand what can be learned about sources that emit 
over different RM ranges. The first two bands are proposed 
for SKAl-mid liPewdnev et al.l2013h . the third band is avail¬ 
able with the Australia Telescope Compact Array (ATCA), 
while the fourth band is available with a number of radio 
telescopes. Each time we will use 1 MHz frequency chan¬ 
nels that have uniform response functions. For each of these 
frequency windows we provide full-width at half maximum 
(FWHM) values of the RM spread function in Table [T] 

Because of the low frequencies and relatively wide fre¬ 
quency channels that we simu late we have to check whether 
we can use the formalism bv iBrentiens fc de Bruv3 (l2005h 
to calculate RM spectra, or the ne w formalism that we pro- 
posed in an accompanying paper llSchnitzeler fc L^l2015l : 
we will refer to this paper as ‘SL15’). In SL15 we show that 
the discrete Fourier transform between wavelength squared 
and RM that is commonly used to calculate RM spectra is 
only approximately correct for most data sets. Exact RM 
spectra can be calculated only if the channel response func¬ 
tion of the data is included in this calculation. The complex 
exponential exp(—2iRM'Ac), where ‘RM'’ is the trial RM 
and Ac = (Af -|- A|) /2 is the average wavelength squared of 
the channel, is exact if frequency channels have a top-hat 
response function in wavelength squared. In SL15 we derive 
how RM spectra can be calculated for frequency channels 
with any type of channel response function, and we pro¬ 
vide expressions for calculating RM spectra if the channel 
response function has a top-hat shape in frequency. Many 
geometries that we simulate emit over a range in RM that 
is much wider than the FWHM of the RM spread function 
(RMSF); to add the contributions by emission at (very) dif¬ 
ferent RMs the wings of the RMSF have to be calculated 
exactly, which requires the formalism that we developed in 
SL15. Therefore we will use in all cases the formalism from 
SL15, in particular equation 9 from that paper and the nor¬ 
malized version of equation 12, to calculate RM spectra from 
the frequency spectra that we simulated. 
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Figure 7. RM spectra for a uniform source with a transverse linear RM gradient from S I2.1.1I which emits over an RM range of 10, 100, 
1000, and 2500 rad (red to blue). The frequency window used in the computation of RM spectra is indicated in the corner of each 
panel. The vertical dotted lines indicate the FWHM of the RM spread function as a crude way for estimating when the RM spectrum 
is resolved. The black spectrum shows the level of instrumental polarization if the source has an intrinsic polarization of 5% and the 
polarization purity of the instrument is -30 dB of Stokes I. The step size in the RM spectra is equal to ten samples per FWHM of the 
RM spread function for the three lowest frequency windows, and twenty samples per FWHM for the window at the highest observing 
frequencies. 


3.1 Instrumental Polarization 


Sources that are severely depolarized require long integra¬ 
tion times to be detected. In addition, polarization leak¬ 
age from the instrument restricts which polarized sources 
can be detected even if the observing time is sufficient. We 
model the contribution of instrumental polarization to the 
RM spectrum as an RMSF centred on RM = 0 rad m~^ 
with a peak amplitude of -30 dB times the total intensity 
signal. We will assume that the source has an instrinsic po¬ 
larization percentage po = 5%, which produces an instru¬ 
mental polarization response with a height of 20 flux den¬ 
sity units. The SKA baseline design specifies a polarization 
purity of -30 dB acro ss the FWHM of the primary beam 
llDewdnev et al.l|2013l ). This performance is better than the 
instrumental polarization level of either the Dominion Ra¬ 
dio Astrophysical Observatory (DRAG) synthesis telescope 
and th e Very Large Array (VLA) prior to its u pgrade to the 
JVLA llTavior et aLll200UlCondon et al.lll998l ). 


3.2 Synthetic RM spectra 

3.2.1 Geometries with large-scale magnetic fields 

In Figs [7] and [8] we show RM spectra that we calculated 
for the uniform source (§ 12.1.11) and the Gaussian source 
(§ 12.1.211 for four frequency windows. As we discussed in 
the previous Section, the frequency spectra of the sources 
that we modelled can be approximated by the spectrum of 
a Gaussian source, and at low and intermediate frequencies 
sometimes by the spectrum of a uniform source. We show 
the FWHM of the RM spread function as a rough indication 
for when the RM spectrum of a source becomes resolved. 

Comparing these two models has the additional advan¬ 
tage that we can investigate how RM synthesis handles the 
sharp edges of the RM spectrum of the uniform source, while 
the RM spectrum of the Gaussian source has smooth edges. 
As SL15 showed, for moderate RMs RM synthesis can be 
approximated w ell by the discrete Fourier tr ansform that 
was proposed bv iBrentiens fc de BruvnI (l2005fl . The rectan¬ 
gular window function then acts as a highpass filter in this 
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Figure 8. RM spectra for Gaussian sources with RM gradients of 10, 50, 125, and 250 rad m~^ FWHM“^ (red to blue), shown in an 
equivalent way to Fig. [7] We show a smaller range in RM along the x-axis because the Gaussian sources that we modelled emit over a 
smaller range in RM than the uniform sources from Fig. [7] 


Fourier transform, and structure on large RM scales will be 
missing from the reconstructed RM spectrum of the uni¬ 
form source, similar to the ‘missing short spacings’ problem 
in radio interferometry. Because the Gaussian source that 
we modelled has a smooth RM spectrum its reconstructed 
RM spectrum suffers much less from this effect. In the low¬ 
est frequency window the uniform source and the RM spread 
function become ragged far from 0 rad m“^, the mean RM 
of the models. This is the result of the discrete frequency 
sampling and the frequency channel response functions that 
we simulated. 

Observations in the lowest frequency band provide the 
most accurate RM measurements, which makes it possible 
to correctly identify sources that emit over RM ranges that 
are only slightly different. Fig. [7] illustrates this for exam¬ 
ple for the two sources that emit over a range in RM out 
to ± 10 rad m~^ and ± 100 rad m“^, which are difficult 
to tell apart from observations in the 1300-3100 MHz win¬ 
dow, but can be identified correctly from observations in 
the two windows at lower frequencies. If sources that emit 
over a wider range in RM have a lower peak polarized flux 
density (equation [5l, then our simulations show that at low 
frequencies only sources that emit over a small range in RM 
can be detected above the instrumental polarization level. 


Uniform sources that emit over a wide range in RM might 
show up only as two peaks above the instrumental polariza¬ 
tion threshold. This effect is less severe at higher frequen¬ 
cies; therefore we recommend including high-frequency data 
(> 1 GHz) to search for sources which emit over a range 
in RM that is larger than several tens of rad m~^. The 5- 
7 GHz band is suitable for identifying sources which emit 
over RM ranges of more than several hundred rad m“^. The 
key role of this frequency band is to determine the amount 
of wavelength-independent depolarization, because Faraday 
rotation is much less severe at high frequencies. 

3.2.2 Geometries with turbulent magnetic fields 

We investigate how RM spectra can help us understand the 
properties of the turbulent magnetic field that we simulated 
in § [U (its strength, coherence length, and the number of 
layers with turbulent cells) in three ways: first, using the 
peak flux density in the RM spectrum, second, the shape 
of the RM spectrum, and third, the difference between the 
RM spectra of partia l coverage mo d els an d our random walk 
models. Previously, iBernet et al.l ||2012|) investigated how 
RM synthesis can be used to investigate the properties of 
turbulent Faraday screens for a small number of turbulent 
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Figure 9. Peak polarized flux density in the RM spectrum versus the largest possible RM for a single turbulent cell for the geometry 
described in § 12.21 for different frequency windows. In each panel the different colours and plot symbols indicate different values for the 
magnetic field strength: 1, 5, 10, 25, and 50 fiG. The vertical dotted line indicates the half-width at half maximum of the RM spread 
function from Table [T] When following a curve that connects points with the same colour and plot symbol the coherence length of the 
magnetic field decreases when going from right to left and the number of sightlines through the Faraday screen increases. We ran 2000 
Monte Carlo simulations to determine the distribution of peak flux densities in the RM spectrum. The plot symbol indicates the median 
value of this distribution, and the error bars enclose the most compact 95% confidence interval. 


cells and a Gaussian pdf(RM). We consider models where 
the number of sightlines through the turbulent foreground 
varies over a much wider range, and we consider different 
probability density functions. 

First, we investigate how the peak normalized flux den¬ 
sity in the RM spectrum depends on the coherence length 
and strength of the magnetic field. Equationexpresses how 
polarized emission from the background source is divided 
into Ni.o.s. polarization vectors, one vector for each line of 
sight through the turbulent foreground medium. These po¬ 
larization vectors all have the same length, equal to P/Ni,o.s., 
but different RMs. The frequency and RM spectra combine 
the behaviour of each of these polarization vectors; the RM 
spectrum consists of Ni,o.b. peaks spread over a range in RM 
between Mayers X RMmax, where Mayers indicates the num¬ 
ber of layers of cells in the turbulent foreground screen, and 
RMmax is defined by equation IIII 

In Figures [5] and [TO] we show the peak polarized flux 
density in the RM spectrum for a number of magnetic field 
strengths, values of RMmax, and two different thicknesses 


of the Faraday-rotating screen. The background source is 
identical to the source that we simulated in §111 it has a 
diameter of 25 parsec and illuminates the Faraday screen 
uniformly. The electron density in the foreground screen is 
10 cm“® and the magnetic field strength is 1, 5, 10, 25, or 50 
nG. We considered a Faraday screen that consists of eight 
layers because its pdf(RM) is very close to Gaussian (Ap- 
pendix|C|. The four panels in these figures correspond to the 
four frequency windows that we used in S 13.2.21 For curves 
of a given colour or given symbols, the only variable is the 
coherence length of the magnetic field in the Faraday screen: 
the coherence length decreases when going from right to left 
in these panel as the number of sightlines Mi.o.s. increases. 
The smallest number of sightlines that we considered is two 
since one sightline produces no depolarization. For magnetic 
fields of 25 and 50 fiG two sightlines produce very large 
RMmax, therefore, instead of starting with two sightlines 
we started with that number of sightlines which produces 
RMmax = 2500 rad . We stop increasing the number of 
sightlines either when the simulation reached RMmax = 10 
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Figure 10. The equivalent of Fig. |9] when there are eight layers of turbulent cells. Note that the maximum RM of a single line of sight 
through the Faraday screen is eight times the maximum RM per cell which we plot along the x-axis of the panels. 


rad m~^ or when the number of sightlines becomes larger 
than 5000. 

To understand the behaviour of the points shown in 
Figs [5] and [TO] consider a single Bturb and let vary. 

When Ni.o.b. is small, RMmax (oc 1/VM.o.s.) is large, and 
since there are not many cells, the polarized flux density 
that is associated with each peak in the RM spectrum 
(= 1000/Afi.o.s. ) is large. The mean separation between two 
peaks in the RM spectrum is equal to 


(ARM) 


2RM,^ax 

iVl.o.s. 


RturbfVla 
^ 3/2 
1.0.s. 


(14) 


When Ni,o.b. increases the length of the polarization vector 
associated with each sightline decreases, and the peaks move 
closer together in the RM spectrum because the RM range 
over which they are distributed (|RM| $1 RMmax) decreases. 
It becomes more likely that some of the peaks in the RM 
spectrum overlap, starting with a partial overlap of the RM 
spread functions. This explains why the peak polarized flux 
density in the RM spectrum does not decrease monotoni- 
cally when W.o.s. is increased. Increasing W.o.s. even further 
moves the peaks closer together, and as a result of over¬ 
lapping RM spread functions the peak flux density in the 
RM spectrum will rise again. The intrinsic position angles 
of the emitted radio waves, xo, are the same in our sim¬ 


ulation, therefore the polarization vectors of the individual 
peaks will align almost perfectly once RMmax becomes much 
smaller than the FWHM of the RM spread function, which 
leads to very high peak polarized flux densities in the RM 
spectra. 

From comparing the four panels in Figs and 1101 and 
matching panels in both Figures, we draw the following con¬ 
clusions. 

(i) If W.o.s. is small (the field coherence length is large 
compared to the size of the background source) and the 
peaks in the RM spectrum do not overlap then models with 
different Bturb and Aiayera produce the same peak polarized 
flux density in the RM spectrum, albeit at different RMmax- 
When Ni,a.B. is small the RM spectrum is not sampled well, 
therefore it is difficult to determine either RMmax or the 
shape of pdf(RM) and from that Mayers- 

(ii) Increasing Bturb whilst keeping Miayers fixed increases 
the RMmax where the peak polarized flux density in Figs [9] 
and [TO] is at its lowest- The minimum in these Figures de¬ 
pends on the mean separation between the peaks in the RM 
spectrum relative to the FWHM of the RM spread func¬ 
tion. To keep the mean separation between peaks in the 
RM spectrum (equation 1 1411 fixed when Rturb is increased, 
Ni.o.b. only has to increase by a small amount. Since RMmax 
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Figure 11. Variances calculated for a Faraday screen with var¬ 
ious number of sightlines passing through them, for Faraday 
screens consisting of 1, 2, and 8 layers of turbulent cells. The 
maximum RM for a single turbulent cell is equal to 50 rad m“^ 
divided by the number of layers, which keeps the maximum total 
RM through the screen fixed and independent of the number of 
layers. Diamonds indicate the median variance, error bars indi¬ 
cate the most compact 95 per cent confidence interval, and the 
horizontal dotted lines indicate the variance that we calculated 
analytically for each of the Faraday screens. 


oc Bturb/VM.o.s. increasing Bturb decreases the RMmax of 
the minimum in FigslOland llOl 

(iii) If Ni.o.s. is large, models with different Bturb pro¬ 
duce the same peak flux density in the RM spectrum at the 
same RMmax. It is then impossible to determine Rturb from 
the measured peak polarized flux density even if RMmax 
is known. However, if M.o.s. is large pdf(RM) is sampled 
well, and it might be possible to determine both RMmax 
and Mayers from the shape of pdf(RM). 

(iv) Sources with a small field coherence length (such that 
Ni.o.s. is large) could be invisible in low-frequency windows 
and only appear in high-frequency windows. This is because 
the peak polarized flux density in the RM spectrum can 
be higher in high-frequency observing windows than in low- 
frequency windows. In low-frequency windows, if A^i.o.s. is 
large the peaks in the RM spectrum have a small ampli¬ 
tude, and the peaks in the RM spectrum do not overlap. At 
higher frequencies the RM spread functions are wider, and 
the peaks start to overlap in the RM spectrum. Because we 
use the same intrinsic position angle xo for all polarization 
vectors, this overlap can increase the peak polarized flux 
density in the RM spectrum. This observation could lead 
to a trade-off between being able to detect a source (high 
peak polarized flux density in the RM spectrum caused by 
overlapping RM spread functions) versus being able to accu¬ 
rately measure the RMs of the individual peaks in the RM 
spectrum (which requires narrow RM spread functions). 

(v) If the RM spectrum is barely resolved and M.o.s. is 
small, as is the case for the 5-7 GHz window for Rturb = 
1 nG, the peak polarized flux densities show complex be¬ 
haviour. 

Second, if pdf(RM) is sampled sufhciently finely it 
should be possible to determine RMmax and the number 
of layers in the Faraday screen from the shape of the RM 


spectrum. We estimate the number of sightlines that is re¬ 
quired to do this using a Monte Carlo simulation, drawing a 
number of sightlines from pdf(RM) for a Faraday screen that 
consists of 1, 2, or 8 layers of turbulent cells, and repeating 
this process 2000 times. The pdf(RM) of a Faraday screen 
which consists of eight layers of cells is almost identical to a 
Gaussian pdf(RM). Fig. 1111 shows the median variance, the 
95 per cent confidence interval of the variances, and the an¬ 
alytically calculated variance for these three different Fara¬ 
day screens. The range in RM that is spanned by pdf(RM) 
is the same for the three Faraday screens, only the shape of 
the distribution is different. As can be expected, if the num¬ 
ber of sightlines is small pdf(RM) is not well sampled, and 
the scatter in the variances for the different Faraday screens 
is large: it is impossible to identify screens of different thick¬ 
ness based on the shape of their RM spectrum. Once there 
are at least several tens of sightlines it becomes possible to 
identify the shape of pdf(RM) and from that the number of 
layers in the Faraday screen. In our simulation we did not 
include the finite width of the RM spread function. Because 
of this finite width the peaks in the RM spectrum from the 
individual sightlines will overlap, which complicates finding 
the shape of pdf(RM) of the Faraday screen. 

Third, and final, we consider how RM spectra of tur¬ 
bulent foregrounds that we modelled in § 12 .2.21 differ from 
RM spectra of partial coverage models. In a partial cover¬ 
age model the Faraday-rotating screen consists of a con¬ 
tribution by a Gaussian pdf(RM) and a contribution by 
the background source that does not undergo Faraday ro¬ 
tation ('equation 1121) . The latter contribution shows up in 
the RM spectrum as a peak at RM = 0 rad m“^, which 
combines with the Gaussian pdf(RM) from sightlines that 
pass through the Faraday screen. The relative heights of the 
Gaussian pdf(RM) and the single peak that is produced by 
sightlines that do not pass through the Faraday screen de¬ 
pends on the fraction of sightlines which pass through the 
Faraday screen. For a source that emits uniformly across its 
surface, this fraction is equal to the fraction of the surface 
of the source that is covered by the Faraday screen. Because 
instrumental leakage can show up as a strong signal at RM 
= 0 rad m“^, telescope leakages should be calibrated well 
to distinguish between the peak in polarized flux density at 
RM = 0 rad m~^ that is produced by the partial coverage 
model and the same peak that is produced by polarization 
leakage in the telescope. 


4 SUMMARY 

We have modelled sources with large-scale and turbulent 
magnetic fields, and predicted frequency spectra for these 
sources between 200 MHz and 10 GHz, where we included 
the finite width of the frequency channels in our analysis. 
In these models Faraday rotation takes place in front of the 
sonrce of the emission. We considered uniform and Gaussian 
sources on the sky with linear transverse RM gradients, and 
cylinders and spheroids with azimuthal magnetic fields. The 
cylinders and spheroids can be inclined with respect to the 
plane of the sky, and for the spheroids we considered differ¬ 
ent axis ratios ranging from 1:1:0.25 (oblate), 1:1:1 (sphere), 
to 1:1:4 (prolate). These source types with large-scale mag¬ 
netic fields do not show net Faraday rotation: this is because 
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integrating over the surface of the source combines sightlines 
with the same polarized flux density but RMs of opposite 
sign. As a result any net RM of the source must be pro¬ 
duced in the foreground, further away from the source. We 
derived the probability density functions (pdfs) of RM for 
turbulent screens with between one and twenty layers of tur¬ 
bulent cells. A single layer of cells has a uniform pdf(RM), 
while a Gaussian pdf(RM) is a good approximation for a 
Faraday screen that consists of more than about four layers 
of turbulent cells. We show how increasing the thickness of 
the Faraday-rotating layer changes pdf(RM) and the prop¬ 
erties of the polarized flux density spectra as a function of 
frequency and RM. 

At high frequencies all the source types which we con¬ 
sidered show a similar drop-off in polarized flux density with 
decreasing frequency. Some source types show secondary 
maxima in their polarized flux density spectra at low and 
intermediate frequencies, which helps with their identifica¬ 
tion. _ 

The model bv iBurnI ll 19661 ') is often used to explain de¬ 
polarization in turbulent foregrounds. It requires that many 
sightlines pass through the Faraday screen, which means 
that the field coherence length is small compared to the ex¬ 
tent of the background source. This large number of sight¬ 
lines in Burn’s model also produces complete depolariza¬ 
tion at low and intermediate frequencies. Partial coverage 
models were developed to explain why some sources are 
not completely depolarized at these frequencies. We show 
that a Monte Carlo model of a small number of sightlines 
passing through a turbulent Faraday screen predicts a drop¬ 
off at high frequencies similar to the Burn depolarization 
model, and a plateau in polarized flux density at low and 
intermediate frequencies. The polarized flux density of this 
plateau can be calculated from a random walk of the po¬ 
larization vectors. Spectra that are produced by the Monte 
Carlo model are very similar in shape to spectra produced by 
partial coverage models; therefore, Monte Carlo models can 
be considered as an alternative to partial coverage models. 

We calculate RM spectra for four frequency windows: 
350-900 and 950-1760 MHz, and 1.3-3.1 and 5.0-7.0 GHz. 
At low frequencies, sources are strongly depolarized and 
might not be detected above the instrumental polarization 
level. However, with sufficient sensitivity RMs can be de¬ 
termined accurately. Similar to the missing short-spacing 
problem in interferometry, if only low-frequency data are 
used the reconstructed RM spectrum of the uniform source 
shows a local minimum at its centre. Such sources can show 
up in RM spectra as two peaks above the instrumental po¬ 
larization level. 

Each sightline through a turbulent Faraday screen pro¬ 
duces a single peak in the RM spectrum; the reconstructed 
RM spectrum is a superposition of these peaks convolved 
with the RM spread function. We analyse the peak height 
and shape of the RM spectra when we change the strength 
and coherence length of the turbulent magnetic field and 
the number of layers with turbulent cells. Decreasing the 
field coherence length increases the number of sightlines and 
therefore the number of peaks in the RM spectrum. At the 
same time the polarized flux of the background source is dis¬ 
tributed over more peaks, and the range in RM over which 
the peaks are distributed is reduced. The combined effect 
of this is that the peak amplitude in the RM spectrum first 


decreases when the field coherence length is decreased, then 
increases when the individual peaks in the RM spectrum 
(convolved with the RM spread function) overlap. When 
more then several tens of sightlines pass through a thin, 
turbulent foreground screen the shape of pdf(RM) can be 
determined reasonably accurately, and with this informa¬ 
tion the number of layers of turbulent cells can be derived. 
Finally, in a partial coverage model the part of a source that 
is not covered by a turbulent screen produces a peak in the 
RM spectrum at 0 rad m~^. Random walk models do not 
show such a peak; therefore RM spectra can tell which of 
the two models applies. 

The modelling framework we provide can be extended 
to any frequency coverage, frequency range, and channel 
width of the observations, and can be used to explore more 
complex source types and combinations of source types. 
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APPENDIX A: FREQUENCY SPECTRA OF A 
CYLINDER WITH AN AZIMUTHAL 
MAGNETIC FIELD 

In this Appendix we model the frequency spectrum of an 
emitting cylinder with a magnetic field that wraps around 
the circular cross-section of the cylinder. Faraday rotation 
occurs in a boundary layer between the inner, emitting cylin¬ 
der and a coaxial outer cylinder. The magnetic field inside 
the emitting cylinder Rem has a constant strength, and 
points along the major axis of that cylinder. The inner cylin¬ 
der has a radius R, while the outer cylinder has a radius R'. 
We define a Cartesian coordinate system {x,y,z) where x 
points towards the observer, y lies along the projection of 
the short axis of the cylinder on the sky, and 2 along the pro¬ 
jection of the long axis. We will use 9 for the angle between 
the plane of the sky and the major axis of the emitting cylin¬ 
der. The observed monochromatic net polarization vector is 
given by equation [2] The E vector of the emitted polar¬ 
ization vector Pern lies along the y axis of the coordinate 
system that we chose. We define the intrinsic position angle 
Xo to be equal to zero degrees if Pem points in the direction 
of y. 

Because in our model the synchrotron-emitting and 
Faraday-rotating layers are not mixed we can solve equa¬ 
tion [2] by first calculating the emitted polarized flux density 
along each line of sight, then calculating the amount of Fara¬ 
day rotation of this emission. To simplify our analysis we will 
not consider sightlines that pass through the polar caps of 
the cylinder. 

The monochromatic volume emissivity for synchrotron 
radiation, depends on the strength of the magnetic 
field inside the cylinder and the inclination of the cylin- 
der relative to the sky 9 as = K |i?em cos 9\°‘^^ i'~°‘ (e.g., 
iRvbicki & LightmanlllQT^ . ‘K’ depends on the mass and 
charge of the synchrotron emitting particles and on the spec¬ 
tral index of the synchrotron emission a {Si, oc The 

flux density of a single line of sight is found by integrating 
the volume emissivity along the line of sight: 

Pem {y, ^) = J (®> V’ P \/l - (y/Rf / |cos 9\ da;, 

where the factor — {y/R)^ corrects for the angle between 
the line of sight and the near surface of the cylinder, and 
1 / |cos0| corrects for the inclination of the cylinder with re¬ 
spect to the plane of the sky. 

In our model the frequency dependence of the emis¬ 
sivity v~°‘ is absorbed into the product of the Stokes I 
spectrum and the wavelength-independent polarization frac¬ 
tion po; by requiring that the source emits 1000 units 
of polarized flux density independent of 9 the entire fac¬ 
tor K |i?em cos / |cos 6| is fixed. Therefore the volume 
emissivity is constant throughout the emitting cylinder, and 
the polarized flux density of a single line of sight depends 
only on the geometry of the emitting cylinder with respect 
to the line of sight: 

(y. „ 1 _ (]lY 

Pol \r) ' 

The amount of Faraday rotation of a single sightline 
can be calculated by integrating the product of the mag¬ 
netic field strength along the line of sight and the local free 
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Figure Al. The behaviour of Stokes Q (with arbitrary units) 
with position along the minor axis of the cylinder, for RMmax = 
1000 rad and frequencies of 3.0 GHz (dark grey) and 1.0 GHz 
(light grey). The Faraday-rotating boundary layer has a thickness 
1/lOth of the radius of the synchrotron-emitting cylinder. This 
profile does not depend on 6 because of the normalization we 
used. 


electron density in the boundary layer between the inner 
and outer cylinder: 


RM(2/,2) = 


^ I 




yfW- 


O.SlnelBI 


= 0 . 81 ne|-B|j/ln 


'R' + ^{R'f-y2 


zda; 


(Al) 


where the boundary layer has an electron density rie and 
is threaded by a magnetic field of amplitude |B|. If the 
source is inclined relative to the plane of the sky then 
the pathlength through the Faraday screen increases as 
1 / |cos0|, which exactly cancels the projection of the mag¬ 
netic field along the line of sight in the Faraday-rotating 
layer (oc |cos0|). Therefore RM(y, 2 ) does not depend on 
the inclination of the source relative to the sky. 

The maximum RM, RMmax, is found at y = R, and 
equation can be expressed in terms of RMmax and the 
geometry of the system as 


RM {y, z) 

RMmax 


In 


R'+V'(R')^-y^ 

R+\/R^-y^ 


R 


In 


f+v/(¥f 


-1 


(A2) 


If R' is only slightly larger than R equation IA2I gives al¬ 
most identical results as when we replace the line-of-sight 
component of the magnetic field, which varies in amplitude 
along the line of sight, with the value halfway through the 
Faraday-rotating boundary layer: 


RM (y, z) 
RMmax 




(A3) 


In Fig. lAll we show the observed values of Stokes Q for 
RMmax = 1000 rad m~^ as a function of y for frequencies of 
1.0 and 3.0 GHz. Our definition for xo = 0° results in the 


emitted polarization vector Pem showing up only in Stokes 
Q, Stokes U is zero. The cylinder that we simulated does not 
show any net Faraday rotation: the symmetry of the system 
is such that points on opposite sides of the major axis of the 
cylinder (at ±i/) emit the same polarized flux density and 
are Faraday rotated in opposite directions, which cancels 
Faraday rotation for the source as a whole. Therefore all net 
polarization vectors for the source as a whole will show up 
only in Stokes Q, independent of the observing frequency. 

To calculate the monochromatic polarization vector for 
the source as a whole we integrated over the minor axis y of 
the cylinder using Romberg’s method. To integrate across 
the finite width of the individual frequency channels, and to 
correctly include changes in the orientation of the monochro¬ 
matic net polarization vector within each frequency channel, 
we divide each channel into subintervals that are narrow 
enough so that we can calculate the net polarization of each 
subinterval using the trapezium rule. We use forty subin¬ 
tervals for each 27t revolution of the position angles that is 
induced by the largest RM of the geometry that we mod¬ 
elled, RMmax. By combining the net polarization vectors of 
the subintervals we calculate the net polarization vector of 
a single frequency channel as a whole. We only include a 
frequency channel in our analysis if the length of the polar¬ 
ization vector of the first subinterval is more than ten times 
the numerical accuracy e defined by equation [B] otherwise 
we stopped calculating the frequency spectrum altogether. 
The first subinterval lies at the low-frequency-end of each 
frequency channel. 


APPENDIX B: FREQUENCY SPECTRA OF AN 
ELLIPSOID WITH AN AZIMUTHAL 
MAGNETIC FIELD 

In this Appendix we model the frequency spectrum of a 
synchrotron-emitting ellipsoid. This geometry consists of 
two nested ellipsoids with a shared coordinate system, shown 
in Fig. [ml The boundary of the inner, synchrotron-emitting 
ellipsoid is described by (x'/A)^ -I- (y'/B)^ + (z' jC'f' = 1, 
and the magnetic field inside the emitting ellipsoid points 
along the z' axis. Faraday rotation occurs in a layer between 
the inner ellipsoid and an outer ellipsoid whose surface is 
given by (x'/Mb)^ + (i/VRb)^ + (zV^'b)^ = 1. The magnetic 
field in the Faraday-rotating layer wraps around the major 
axis {z' axis) of the inner ellipsoid and is parallel to the 
surface of the inner ellipsoid. The thickness of the Faraday- 
rotating boundary layer varies because the Faraday-rotating 
layer is nested between two ellipsoids; by comparison, in our 
model of the emitting cylinder (Appendix]^ the thickness 
of the Faraday-rotating boundary layer was the same every¬ 
where. 

In the [x' ,y' ,z') coordinate system of the ellipsoid the 
line of sight towards the observer, indicated by the unit vec¬ 
tor nios, is specified by the angles ij) a-nd i/’, where 

nios = {cos ?/), sin ?/) cos (j), sin 'ip sin (p} , 

We choose the two mutually perpendicular unit vectors 
and ?/p in the plane of the sky as the coordinate basis, where 

Xp = {—sin-i/), cos <()cos i/’, sin <()cos i/’, } , 

y? = (0, — sin (p, cos (p} . 
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Figure Bl. Geometry relating the coordinate system of the el¬ 
lipsoid with surface S and magnetic field B to the line of sight 
‘los’ towards the observer ‘obs’ .The orthonormal vectors aip and 
j/p lie in the plane of the sky. 


In the main text we will use (j) = 90°. In that case 'tp = d, the 
inclination of the major axis of the ellipsoid with respect to 
the plane of the sky. 

To calculate the monochromatic net polarization vector 
P {v), which is integrated over the surface of the source, we 
divide the surface of the source into a very fine grid of sight¬ 
lines and calculate the polarization vector of each of these 
sightlines numerically. Similar to the polarized emission from 
the cylinder that we considered in Appendix we assume 
that both the inclination angle of the magnetic field and the 
volume emissivity are uniform inside the emitting ellip¬ 
soid. The inclination angle of the magnetic field can then 
be absorbed into the normalization of the monochromatic 
net polarization vector, and the emitted polarized flux den¬ 
sity is simply proportional to the length of the line of sight 
through the emitting ellipsoid. We calculate the polarized 
flux density of each sightline analytically; the proportion¬ 
ality constant is fixed by our requirement that the source 
emits 1000 units of polarized flux density. 

To calculate the amount of Faraday rotation of each 
sightline we integrate equation [T] numerically through the 
Faraday-rotating boundary layer using the trapezium rule, 
with 40 sampling points per sightline. Then we calculate the 
monochromatic net polarization vector P (v) by adding the 
polarization vectors for a square grid of 1200 x 1200 sightlines 
across the surface of the source projected onto the sky. For a 
few of our models we varied the number of sampling points 
through the Faraday rotating layer (20 or 40 points) or the 
density of the grid of sightlines across the surface of the 
source (800x800 or 1200x1200 sightlines). By comparing 
the results from these models we found that the fractional 
accuracy in the length of the polarization vector is better 
than 1% of the emitted 1000 units of polarized flux density. 
Close to the major axis of the source, where most of the 
polarized flux is generated, each 27t wrap of the position 
angles across the surface of the source is sampled with at 
least 2 sightlines. 

To integrate P (v) across the finite width of the fre¬ 
quency channels we calculate P (v) at forty points for each 
27t wrap in position angle that is induced by the sightline 


with the largest RM, and we integrate across each frequency 
channel using the trapezium rule. 


APPENDIX C: THE PROBABILITY DENSITY 
FUNCTION OF RM OF A TURBULENT 
FOREGROUND MEDIUM 

If there is only one layer of turbulent cells in front of the 
emission region then the RMs from this layer follow a uni¬ 
form distribution. This can be shown as follows. We will use 
‘pdf’ as shorthand for the probability density function, and 
i indicates the angle between the line of sight and the mag¬ 
netic field. The observed RM = RMmaxCosi, where RMmax 
is defined in equation Because probability mass is con¬ 
served, 

odfiRMl = = l/2\smi\ ^ 1 

|dRM/di| RMmax |sini| 2RMmax ’ 

therefore pdf(RM) is uniform between ± RMmax. pdf (i) is 
equal to the fraction of the surface on the unit sphere be¬ 
tween i and i + di. 

For two layers of turbulent cells the pdf for measuring 
RMi ,2 = RMi + RM 2 , pdf (RMi, 2 ), can be found as the 
convolution pdf (RM) * pdf (RM), 

/ OO 

pdf (RMi,2 - RM) pdf (RM) dRM, (Cl) 

-00 

which is triangular. For three layers, pdf (RMi, 2 , 3 ) = 

{pdf(RM)*pdf(RM)} * pdf(RM) = pdf(RMi,2) * 
pdf (RM). This way pdf (RMi,.,,n) can be calculated for any 
number N of layers with turbulent cells. 

Fig. m shows pdfs for between one and twenty lay¬ 
ers of turbulent, Faraday-rotating cells and their matching 
normal distributions, while Fig. IC2I shows the difference in 
probability density between the normal distributions and the 
pdf(RM). The standard deviation a of each normal distribu¬ 
tion follows from the variance of probability density function 
that we calculated analytically. The maximum difference be¬ 
tween the normal distribution and pdf(RM) decreases from 
0.0380 for 2 layers of Faraday-rotating cells to 0.0071 for six 
layers, 0.0033 for ten layers, and 0.0012 for twenty layers of 
cells (these are the values for the pdfs themselves, not for 
pdf/max(pdf) that we show in Fig. IcB . 
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Figure Cl. Probability density functions for one, two, four, six, 
ten, and twenty layers of turbulent cells (grey lines) and the nor¬ 
mal distributions found by calculating the standard deviation of 
the probability density functions (dashed lines; see the text for 
details). 



Figure C2. Difference in probability density between scaled nor¬ 
mal distributions from Fig. ICll and the pdf(RM), for different 
numbers of layers with turbulent Faraday-rotating cells. 








